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I.  INTRODUCTION 

At  the  last  symposlun  Schaaf  and  Maslach  (1963)  presented  experimental 
cylinder  drag  results  which  were  obtained  at  sufficiently  high  Knudsen  numbers 
and  with  small  enough  scatter  in  the  data  to  make  a  meaningful  comparison  with 
the  theories  then  available  (Lunc  and  Lubonski,  1956;  Baker  and  Charwat,  1958; 
Willis,  1959)  for  high  Knudsen  number  flows.  A  large  discrepancy  between  theory 
and  experiment  was  noted  at  that  time.  As  theoretical  results  were  not  then 
available  for  a  cylinder  the  comparison  was  basec  on  results  for  a  two-dimensional 
strip  normal  to  the  flow.  However,  it  was  felt  that  this  geoinetric  effect  would 
not  significantly  change  the  qualitative  nature  of  the  comparison. 

Since  the  last  meeting  the  following  steps  have  been  taken  in  an  attempt 
to  reconcile  the  results  of  theory  and  experiment.  Firstly,  new  experimental 
results  have  been  obtained  (Section  J.I).  Secondly,  theoretical  calculations 
have  been  performed  for  the  cylinder  (Taub,  1964).  Thirdly,  a  method  different 
from  that  of  Schaaf  and  Maslach  (1963)  has  been  used  to  relate  the  Knudsen 
number  to  the  parameter  which  arises  in  those  theoretical  calculations  where  a 
statistical  model  is  used  to  represent  the  intennolecular  collision  process. 

This  is  discussed  in  Section  III. 

II.  NEW  EXPERIMENTAL  RESULTS 

Since  1962  a  continuing  program  of  drag  force  measurement  in  the  near  free 
molecule  and  free  molecule  regimes  has  been  carried  out  at  Berkeley  extending 
the  range  of  flow  conditions  to  higher  Mach  numbers  and  utilizing  both  cylin¬ 
drical  and  flat  strip  models.  Essentially,  the  same  experimental  methods  have 
been  used  as  were  previously  described  (Maslach  and  Schaaf,  1963)  with  an 
Improved  balance  mechanism.  Fig.  1,  and  the  use  of  free  jet  testing  techniques 
to  achieve  higher  Mach  numbers  at  lower  densities.  The  extensive  precautions 


taken  to  determine  precisely  the  cliaracteristics  of  the  flow  fields  issuing 
from  these  simple  convergent  nozzles  are  the  subject  of  a  separate  paper 
(Ashkenas  and  Sherman,  1964). 

Cylinder  drag  data  were  obtained  at  Mach  numbers  of  approximately  6, 

8,  and  10  with  a  range  of  Knudsen  numbers  from  approximately  1  to  32  (Tang, 
1964).  Flat  strip  drag  data  for  the  model  normal  to  the  flcrw  were  obtained 
at  approximately  the  same  Hach  numbers  with  a  similar  range  of  Knudsen 
numbers  (Ko,  1964). 

All  results  approach  diffuse  free  molecule  flow  limiting  conditions  for 
high  Knudsen  numbers;  typical  data  are  shown  in  Figs.  2  and  3.  PrevioiiS 
results  at  Mach  approximately  5.92  using  uniform  fiv-w  field  testing  tech¬ 
niques  were  essentially  duplicated,  the  maximum  discrepancy  being  2.871, 
utilizing  the  free  jet  testing  methods,  Fig.  4. 


III.  RELATION  OF  KNUDSEK  NUMBER  TO  COLLISIC^  PARAMETERS 

All  calculations  to  date  for  external  aerodynamic  problems  have  used 
either  perfectly  elastic  spheres  to  represent  the  collision  process  (Lunc 
and  Lubonski,  1956;  Baker  and  Cbarwat,  1958)  or  replaced  the  Maxwell -Boltzmann 
collision  operator  with  a  statistical  type  model  (Willis,  1959,  1960;  Taub, 
1964;  Rose,  1964). 

For  hard  sphere  molecules  the  mean  free  path  is  given  unambiguously  by 

/  2 

X  =  (  V  2  n  n  a  )  (1) 

and  the  (Chapman-Enskog)  viscosity  is  given  by 
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Two  statistical  isodels  have  been  used.  Willis  (1959  and  1960)  and 
Taub  (1964)  used  the  ’‘modified  Krook"  model.  This  model  was  specifically 
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designs;?  for  nearly  free  molecular  cond:\totiS.  For  convenience  a  brief 
description  is  given  in  the  Appendix.  Rose  (1964)  uses  the  origi'ial  Krook 
model  (lihatnager,  et .£l . ,  1954) 

Nearly  free  molecular  analyses  using  either  of  these  models  lecclts 

in  terms  of  a  "natursl"  small  parameter  a  =  6n^(D/2)  (m,/2kT^)^'^^ .  In  this 
relation  D  is  a  typical  body  dimension,  n^  the  number  density  at  infinity. 
6  is  a  constant  appearing  in  the  model  (such  that  is  the  relaxation 

time  for  the  translational  degi'ees  of  freedom),  and  is  the  temperature 

characterizing  tlte  Maxwell  distribution  which  the  molecules  reflected  from 
the  body  are  assumed  to  possess. 

In  comparing  theory  with  experiment  we  must  relate  Q  to  the  free  stream 
Rnudsen  number,  Schaaf  and  Maslach  (1963)  used  the  following  relation: 

a  =  0.6  S,/Kn  (3) 

r 

2  2 

where  S  *  ta  /2kT  and  Kn^  «  X This  relation  is  empirical  and  was 
obtained  by  a  comparison  of  high  Hach  number  strip  drag  results  using  the 
modified  Kxook  model  and  hard  sphere  moieculss. 

A  very  different  result  is  obtained  if  we  choose  to  define  mean  free 
path  fr«a  the  viscosity  formula  of  Eq.  (2).  This  formula  is  only  valid  for 
near  equilibrium  conditions,  i.e.,  far  from,  the  body.  For  the  simpK  Krook 
njodel  we  have  the  result 


and  hence 

a  «  (Jt^/74)(T  /T,  O) 

C7  O  00 

There  will  be  a  large  difference  between  the  results  of  Eqs.  (3)  and  (5)  if 
the  speed  ratio,  S  ,  is  not  close  to  unity.  For  the  modified  Krook  model 
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the  Dsear.ing  of  the  vi3cos5.ty  is  not  so  clear  as  the  -sodel  is  specifically 

designed  for  nonequil  ibriursi  'ituations  However,  we  propose  to  use  the 

relation  in  Eq.  (5)  for  this  case  also 

There  are,  of  course,  several  other  ways  in  which  a  could  be  related 

to  Kn^^,  Both  the  physical  arguments  underlying  the  first  collision  methods 

and  order  of  magnitude  analyses  using  the  first  iterate  method  with  the 

model  equations  suggest  that  the  main  •source  of  correction  to  the  free 

molecular  drag,  etc.,  comes  from  collisions  between  molecules  leaving  the 

body  and  free  stream  molecules  A  relation  between  a  and  Kn  might 

00 

t'lcn  be  developed  by  studying  the  rate  of  such  collisions.  This  will  not 
be  act.;.'t<pted  here. 


IV .  CCHFAKIgiftS’  BETWEEN  THEORY  AND  EXPERIMENT 

Unless  specif ically  stated  otherwise  we  use  Eq.  (5)  to  relate  Ci  and 

Kn  for  those  theoretical  results  that  use  the  statistical  models. 

00 

A.  Cylinder  and  ?t^o-D  ?' mens  tonal  Strip  Drag 

The  theoretical  results,  based  on  the  modified  .<rook  model,  are  of 

the  form 


r 
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The  method  used  is  that  of  integral  it  nation  (Willis,  1959),  The  function 
K  is  given  in  Table  I  for  typical  values  of  S  and  S,  used  in  the  drag 

00  O 

expsrin>ents.  It  should  be  noted  that  ter;us  of  order  a  have  been  neglected 
compared  to  those  of  erder  a  in  a.  This  r*a5:urally  limits  the  range  of 
validity  of  the  theory  significantly.  Typical  cosparisons  of  the  data  with 
the  theoretical  results  are  shown  in  Figs.  3  and  5.  For  the  cylinder  we  show 
the  effect  of  the  two  interpretations  of  the  C2-Kn  relation.  It  is  obvious 

CD 

that  the  agreement  's  far  superior  using  the  mean  free  p>.ath  Xased  on  viscos¬ 
ity,  i.e.,  Eq.  (5). 
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Cylinder  Equilibrium  Temperature 

The  .cylinder  equilibrium  censperature  can  be  calculated  from  results 
given  by  Tf.ub  (1964)  for  the  heat  transfer  corrections.  The  calculation 
must  be  performed  using  the  value  of  6.^^^  uniquely  determined  by  S^, 
which  corresponds  to  zero  heat  transfer  for  the  free  molecular  flow.  Using 
the  modified  Krook  model  we  obtain  the  following  result  for  the  normalized 
equilibrium  temperature 


* 
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T  -  T  (Kn  =0) 

-eg _ 91.  “ 


T  (Kn  «« 
eq  o 
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T  (Kn  »0) 
eq  ® 


1  +  (a  in  a)G(S  ,S.  ) 

<x>  D 


(7) 


Assuming  a  monatomic  gas,  a  recovery  factor  of  unity  in  continuum  condi¬ 
tions,  and  supersonic  flow,  we  have  approximately 

T  (Kn  «0)/T  ^[Kn  ««»]  =  (1  +  0.4  S  ^)(T Jt  „[Kn=«])  (3) 

Using  this  result  and  the  standard  formula  for  (T jg/T^q[Kn»«j)  for  a 
cylinder  (Schaaf,  1963)  we  obtain 


s  - 

C9 

2, 

Sb- 

1.12, 

G  »  0,190 

s  ** 

Sb" 

1.31, 

G  “  0.304 

s  ® 

03 

5, 

1.35, 

G  =  0.340 

These  results  are  plotted  in  Fig.  6  using  Eq.  (5)  to  relate  a  to  Kn^. 
The  data  are  taken  from  Dewey  (1961).  Also  shown  is  the  result  given  by 
Schaaf  and  Haslach  (1963  ),  namely 


* 
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,  0  <  86  , *  „  » 
I  ••  - -  (in  Kn  ) 


(10) 


which  was  obtained  from  results  in  (^iilis, 1959)  using  hard  sphere  molecules 
with  S  approaching  infinity.  The  modified  Krook  results  fall  within  the 
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scaccer  bounds  of  f he  data  for  Kn  >5  While  this  agreement  is  scarcely 

CC 

good,  it  is  far  better  than  the  gross  discrepancy  between  the  result  of 
Eq.  (10)  and  the  data 

The  value  of  G  is  quite  sensitive  to  the  assumption  regarding 
the  recovery  factor  under  continuum  conditions  and  may  be  expected  to  vary 
significantly  for  diatomic  gases.  There  is  also  some  question  regarding 
the  Prandtl  number  implied  by  the  model. 

C •  Sphere  Drag 

Sphere  drag  has  been  predicted  by  first  collision  (Baker  and 
Charwat,  1958)  and  first  iterate  (Willis,  1960)  methods.  Recently  Rose 
(1964)  derived  an  expression  for  the  drag  by  a  much  different  analytical 
method  involving  the  use  of  the  linearized  Krook  equation  and  Fourier 
transform  techniques.  All  the  above  results  were  for  large  free  stream 
Mach  number  and  S,  >  2.  Rearranging  all  the  results  in  terms  of  Kn^, 


we  find 


S  ■  Sfm  =  -  (0.24  +  1  06)/to 


(Baker  and  Charwat)  (11a) 


(0,165  S^  +  1.44  -  1.13/S^)/(S  K.n  )  (Willis)  (lib) 

Q  0  CC  00 


(0.33  S,  -  0  12)/S  Kn 

b  CO  c 


(Rose)  (11c) 


Equation  (lib)  is  obtained  by  fitting  numerical  results  for  >  2. 

We  have  compared  the  theories  with  typical  data  given  by  Kins  low 
and  Potter  (1963)  in  Fig.  7.  (We  have  plotted  mean  values  when  mote  than 
one  measurement  was  taade  at  the  same  nominal  conditions.)  The  data  were 
obtained  for  S  =  8.8  and  S,  >  4,4  so  all  the  theories  should  be  applicable. 

CO  0 

Formula  (11a)  is  slightly  different  from  that  given  by  Schaaf  and  Maslach 
(1963)  due  to  a  difference  in  expressing  one  of  Baker  and  Charwat  s  par¬ 
ameters  in  terms  of  Kn  and  S,  ■  Equation  (Ilb)  corrects  a  typographical 

OJ  Q 

error  in  a  corresponding  equation  in  (\tillis.  1960). 
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It  can  be  seen  that  the  results  based  on  hard  sphere  molecules  (Eq 
11a)  are  again  in  serious  disagreement  with  the  data  The  results  for  both 
models,  however,  are  in  quite  good  quantitative  agreement  with  tne  data. 

These  remarks  hold  true  for  all  of  Klnslow  and  Potter's  data  with  =  6.25 
(Fig.  7),  5.85,  5.0,  and  4.4.  The  relative  difference  between  the  results 
of  Eqs.  (lib)  and  (lie)  increases  as  decreases,  but  values  of 

significantly  lower  than  4  4  will  be  needed  to  discriminate  between  the 
results. 

V.  CONCLUSIONS 

1)  New  cylinder  and  flat  strip  drag  data  for  models  normal  to  the 
flow  have  been  obtained  utilizing  free  jet  testing  techniques.  For  similar 
Mach  numbers  the  data  essentially  duplicate  previous  work  using  uniform 
flow  field  testing  techniques. 

2)  For  the  drag  of  a  cylinder  or  two-dimensional  strip  the  data 
have  been  obtained  for  sufficiently  high  Knudsen  number  and  with  sufficiently 
low  scatter  to  provide  a  test  of  the  theoretical  prediction.?. 

3)  The  theoretical  results  obtained  from  the  first  iterate  method; 
using  the  taodified  Krook  model  to  represent  the  collision  process,  and 
determining  the  mean  free  path  from  viscosity  in  the  gas  far  from  the  body, 
are  in  reasonable  agreement  with  experimental  data  for  cylindfr.  two- 
dimensional  strip,  and  sphere  drag.  For  the  cylinder  equilibrium  temper¬ 
ature  the  agreement  is  only  fair,  but  at  least  (T  -I  [Kn  =<»])  has  the 
correct  order  of  magnitude  for  Kn^  >  5. 

4)  Results  obtained  using  the  linearized  Krook  model  (Rose,  1964) 
are  also  in  reasonable  agreement  with  the  sphere  drag  data  of  Kinslow 
and  Potter  (1963). 

I 
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5.  Theoretical  results  using  hard  sphere  molecules  to  represent  the 
collision  process  do  not  agree  with  the  above  cited  experimental  data. 

6.  A  more  definitive  comparison  could  be  obtained  if  the  theoretical 
calculations  for  the  two-dimensional  bodies  were  extended  to  include  the 
terms  of  order  (l/Kn  )  as  well  as  (l/Kn  )  in  (Kn  ),  and  if  the  sphere 

00  CO  00 

drag  results  could  be  obtained  at  higher  Knudsen  numbers  and  for  a  wider 
range  of  body  to  free  stream  temperature  ratios. 


mSLE  I 


values  of  H(S  ,S,  )  (Bq.  6) 
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APPENDIX 

THE  MODIFIED  KROOK  MODEL 


The  simple  Krook  model  (Bhatnager  aj^.  .  I95A)  considers  all 

collisions  at  a  point  statistically  and  makes  no  distinctions  between 
the  various  types  of  collision?  that  occur.  This  appears,  intuitively  , 
to  be  too  simple  a  representati<  for  a  nearly  free  molecular  flow  where 
the  distribution  function  has  very  different  properties  depending  on 
whether  or  not  the  molecular  velocity  is  within  the  solid  angle  subtended 
by  the  body  and  directed  away  from  the  body.  As  a  step  toward  the 
Kaxweli-Boitjnsann  type  of  collision  operator,  where  all  possible  types 
of  collisions  are  considered,  we  propose  the  following  model,  specifi¬ 
cally  designed  for  aonequi librium  conditions.  The  molecules  at  any 
point  r  are  divided  into  two  classes.  Those  whose  molecular  velocities 
(I)  lie  in  the  outward  drawn  solid  angle  subtended  by  the  body  are 
called  class  b  and  all  other  molecules  are  called  class  c.  The  colli¬ 
sion  term  becomes 
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' COillSl 


ions 


=  -  6f,  , +n  w  ) 

b'  o  bb  c  be' 


2  2 
4-  5(ii  w.  .'S’.  . +2n,  n  w,  +n  w  f  ) 
b  bb  ob  0  c  be  be  c  cc  cc 


(Al) 


wtere  6  is  a  constant, 
and 


bb’  be  cb' 


w  are  functions  of  r  only 
cc  — 


exp(.-[m/2kT  JL^-u  ]^)  (A2) 

(The  equation  for  f^  is  given  by  exchanging  b  and  c  in  Eq.  Al). 

The  other  parameters  are  given  by  net  conservation  considerations  as 


FIG.  I  MODEL,  BALANCE  AND  NOZZLE  IN  WIND  TUNNEL 
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FIG.  6  COMPARISON  OF  THEORY  AND  EXPERIMENT, 
CYLINDER  EQUILIBRIUM  TEMPERATURE 
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